Realizable Hamiltonians for universal adiabatic quantum computers 
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It has been established that local lattice spin Hamiltonians can be used for universal adiabatic 
quantum computation. However, the 2-local model Hamiltonians used in these proofs are general 
and hence do not limit the types of interactions required between spins. To address this concern, 
the present paper provides two simple model Hamiltonians that are of practical interest to experi- 
mentalists working towards the realization of a universal adiabatic quantum computer. The model 
Hamiltonians presented are the simplest known quantum-Merlin- Arthur-complete (QMA-complete) 
2-local Hamiltonians. The 2-local Ising model with 1-local transverse field which has been realized 
using an array of technologies, is perhaps the simplest quantum spin model but is unlikely to be 
universal for adiabatic quantum computation. We demonstrate that this model can be rendered 
universal and QMA-complete by adding a tunable 2-local transverse a x a x coupling. We also show 
the universality and QMA-completeness of spin models with only 1-local a z and a x fields and 2-local 
a z a x interactions. 
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What are the minimal physical resources required for 
universal quantum computation? This question is of in- 
terest in understanding the connections between phys- 
ical and computational complexity, and for any practi- 
cal implementation of quantum computation. In 1982, 
Barahona [1] showed that finding the ground state of the 
random field Ising model is NP-hard. Such observations 
fostered approaches to solving problems based on classi- 
cal [2] and later quantum annealing [3] . The idea of using 
the ground state properties of a quantum system for com- 
putation found its full expression in the adiabatic model 
of quantum computation [4] . This model works by evolv- 
ing a system from the accessible ground state of an initial 
Hamiltonian Hi to the ground state of a final Hamilto- 
nian Hi, which encodes a problem's solution. The evo- 
lution takes place over parameters s G [0, 1] as H(s) = 
(l — s)Hi + sH{, where s changes slowly enough that tran- 
sitions out of the ground state are suppressed [5]. The 
simplest adiabatic algorithms can be realized by adding 
non-commuting transverse field terms to the Ising Hamil- 
tonian: V //,^ + V A.rrf + J^afaJ, (c.f. [6]). 
However, it is unlikely that the Ising model with trans- 
verse field can be used to construct a universal adiabatic 
quantum computer [7]. 



What then are the simplest Hamiltonians that al- 
low universal adiabatic quantum computation? For this 
we turn to the complexity class quantum-Merlin- Arthur 
(QMA), the quantum analog of NP, and consider the 
QMA-complete problem k-LOCAL Hamiltonian [8]. One 
solves k-LOCAL Hamiltonian by determining if there ex- 
ists an eigenstate with energy above a given value or be- 
low another — with a promise that one of these situations 
is the case — when the system has at most k-local inter- 
actions. A Yes instance is shown by providing a witness 



eigenstate with energy below the lowest promised value. 

The problem 5-local Hamiltonian was shown to be 
QMA-complete by Kitaev [8]. To accomplish this, Ki- 
taev modified the autonomous quantum computer pro- 
posed by Feynman [9]. This modification later inspired 
a proof of the polynomial equivalence between quantum 
circuits and adiabatic evolutions by Aharonov et al. [10] 
(see also [11, 12]). Kempe, Kitaev and Regev subse- 
quently proved QMA-completeness of 2-local Hamil- 
tonian [14] . Oliveira and Terhal then showed that uni- 
versality remains even when the 2-local Hamiltonians act 
on particles in a subgraph of the 2D square lattice [15]. 
Any QMA-complete Hamiltonian may realize universal 
adiabatic quantum computation, and so these results are 
also of interest for the implementation of quantum com- 
putation. 

Since 1-local Hamiltonian is efficiently solvable, 
an open question is to determine which combinations of 
2-local interactions allow one to build QMA-complete 
Hamiltonians. Furthermore, the problem of finding 
the minimum set of interactions required to build a 
universal adiabatic quantum computer is of practical, as 
well as theoretical, interest: every type of 2-local inter- 
action requires a separate type of physical interaction. 
To address this question we prove the following theorems: 

Theorem 1. The problem 2-local ZZXX Hamil- 
tonian is QMA-complete, with the ZZXX Hamiltonian 
given as: 



tfzzxx = J2 h ^ +J2 A ^ + (!) 

i i 
»J i,3 
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Theorem 2. The problem 2-LOCAL ZX Hamiltonian 
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is Q MA -complete, with the ZX Hamiltonian given as 

Hzx = £>o* I X A << • (2) 

i i 

+ X • / -.- < < + Y, A -,- ■ 

i<j i<j 

a. Structure In the present paper we briefly review 
the standard circuit to adiabatic construction to show 
that 2-local Hamiltonian is QMA-complete when re- 
stricted to real-valued Hamiltonians. We then show how 
to approximate the ground states of such 2-local real 
Hamiltonians by the ZX and ZZXX Hamiltonians. We 
conclude this work by providing references confirming our 
claim that the Hamiltonians in Eq. (1) and (2) are highly 
relevant to experimentalists attempting to build a univer- 
sal adiabatic quantum computer. 



I. THE PROBLEM 

The translation from quantum circuits to adiabatic 
evolutions began when Kitaev [8] replaced the time- 
dependence of gate model quantum algorithms with spa- 
tial degrees of freedom using the non-degenerate ground 
state of a positive semidefinite Hamiltonian: 

= ff|Vhist) = (3) 

(Hi n + -ffclock + "clockinit 

To describe this, let T be the number of gates in the 
quantum circuit with gate sequence Ut • ••U^Ui and let 
n be the number of logical qubits acted on by the circuit. 
Denote the circuit's classical input by \x) and its output 
by iV'out}- The history state representing the circuit's 
entire time evolution is: 
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|a?}<8> \0)® T + U 1 \x)®\l)\Q) 



®T-1 



+ U 2 U x \x)® |11)|0)® T - 2 
+ ... 

+ U T ---U 2 U 1 \x) <g> |1)® T 



(4) 



where we have indexed distinct time steps by a T qubit 
unary clock. In the following, tensor product symbols 
separate operators acting on logical qubits (left) and 
clock qubits (right). 

H m acts on all n logical qubits and the first clock 
qubit. By annihilating time-zero clock states coupled 
with classical input x, H ln ensures that valid input state 
(| a;) <8> 1 0. . .0) ) is in the low energy eigenspace: 



H m = 5^(l-|a*>fci|)®|0><0|i (5) 
=i 

T) x>-(-i) xi o®(i+o- 



-ffciock is an operator on clock qubits ensuring that 
valid unary clock states |00...0), |10..0), 1 110..0) etc., span 
the low energy eigenspace: 



T-l 



iW=]T|01)<01| (M+1) (6) 
t=l 

T-l 

(T - 1)1 + of - <4 - £>"<7ft+l) 



t=l 



where the superscript (t, t + 1) indicates the clock qubits 
acted on by the projection. This Hamiltonian has a sim- 
ple physical interpretation as a line of ferromagnetically 
coupled spins with twisted boundary conditions, so that 
the ground state is spanned by all states with a sin- 
gle domain wall. The term i/ciockint applies a penalty 
|l)(l|t=i to the first qubit to ensure that the clock is in 
state |0)' 8,T - at time * = 0. 

Hp m p acts both on logical and clock qubits. It en- 
sures that the ground state is the history state corre- 
sponding to the given circuit. H prop is a sum of T 

terms, H prop = Ylt=i Hprop,t-> where each term checks 
that the propagation from time t — 1 to t is correct. For 
2 < t < T — 1, i/ pr0 p,t is defined as: 

ffprop.t = f t®\t-l}(t-l\-U t ®\t}(t-l\ 

®\t-l)(t\ + t<Z>\t)(t\, (7) 

where operators \t)(t - 1| = |110)(100| (t— i,t,t+i) stc, act 
on clock qubits t — 1, t, and t+1 and where the operator 
Ut is the t th gate in the circuit. For the boundary cases 
(t = 1,T), one writes H proPtt by omitting a clock qubit 
(t — 1 and t + 1 respectively). 

We have now explained all the terms in the Hamilto- 
nian from Eq. (3) — a key building block used to prove the 
QMA-completeness of 5-local Hamiltonian [8]. The 
construction reviewed in the present section was also 
used in a proof of the polynomial equivalence between 
quantum circuits and adiabatic evolutions [10]. Which 
physical systems can implement the Hamiltonian model 
of computation from Eq. (3)? Ideally, we wish to find a 
simple Hamiltonian that is in principle realizable using 
current, or near- future technology. The ground states 
of many physical systems are real- valued, such as the 
ground states of the Hamiltonians from Eq. (1) and (2). 
So a logical first step in our program is to show the QMA- 
completeness of general real-valued local Hamiltonians. 



The QMA-completeness of real-valued 
Hamiltonians 



Bernstein and Vazirani showed that arbitrary quan- 
tum circuits may be represented using real-valued gates 
operating on real- valued wave functions [17]. Using this 
idea, one can show that 5-local real Hamiltonian 
is already QMA-complete — leaving the proofs in [8] oth- 
erwise intact and changing only the gates used in the 
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circuits. H m from Eq. (5) and -ff c iock from Eq. (6) are al- 
ready real-valued and at most 2-local. Now consider the 
terms in H prop from Eq. (7) for the case of self-inverse 
elementary gates Ut = U}: 

#prop, f = i(l-<7f t _ 1) )(l + af f+1) ) (8) 

For the boundary cases (t = 1,T), define: 
H plopA = i(l + o5)-^®i(af +afoJ) (9) 

H pIO p,T = ~ a (T-l)) ~ U T ® ^i a T ~ a (T-l) a T)- 

The terms from Eq. (8) and (9) acting on the clock space 
are already real-valued and at most 3-local. As an ex- 
plicit example of the gates U t , let us define a universal 
real-valued and self-inverse 2-qubit gate: 

RiM = i(l + <7f) + i(l-<7?)®(sin(#7?+ costal). 

The gate sequence R i j(<p)R i j(ir/2) recovers the universal 
gate from [18]. This is a continuous set of elementary 
gates parameterized by the angle <f>. Discrete sets of self 
inverse gates which are universal are also readily con- 
structed. For example, Shi showed that a set comprising 
the C-NOT plus any one-qubit gate whose square does 
not preserve the computational basis is universal [13] . We 
immediately see that a universal set of self-inverse gates 
cannot contain only the C-NOT and a single one-qubit 
gate. However, the set {C-NOT, X, cos ipX + s'mipZ} is 
universal for any single value of ip which is not a multiple 

of 7r/4. 

A reduction from 5-LOCAL to 2-LOCAL Hamiltonian 
was accomplished by the use of gadgets that reduced 3- 
local Hamiltonian terms to 2-local terms [14] . From the 
results in [14] (see also [15]) and the QMA-complctcncss 

of 5-LOCAL REAL HAMILTONIAN, it now follows that 2- 

LOCAL real Hamiltonian is QMA-complete and uni- 
versal for adiabatic quantum computation. We note that 
the real product af®a^ or tensor powers thereof, are not 
necessary in any part of our construction, and so Hamil- 
tonians composed of the following pairwise products of 
real-valued Pauli matrices are QMA-complete and uni- 
versal for adiabatic quantum computation [29]: 

{1, 1 ® a x ,t ® a z ,a x ® 1, (10) 
a z ® 1, a x ® ct z , a z ® a x ,a x ® a x ,a z ® a z }. 

To prove our Theorems (1) and (2), we will next show 
that one can approximate all the terms from Eq. (10) 
using either the ZX or ZZXX Hamiltonians the Hamil- 
tonians from Eq. (1) and (2) respectively. We do this 
using perturbation theory [14, 15] to construct gadget 
Hamiltonians that approximate the operators of a x and 
a i °f w hh terms from the ZZXX Hamiltonian as well as 
the operators a z a z and erf a x with terms from the ZX 
Hamiltonian. 



B. The ZZXX gadget 

We use the ZZXX Hamiltonian from Eq. (1) to con- 
struct the interaction a z a x from a x a x and a z a z intcrac- 
tions. Let i? c ff = onja z a x ® |0)(0|/-, where qubit k is an 
ancillary qubit and define the penalty Hamiltonian H p 
and corresponding Green's function G(z) as follows: 

H p = 5\l){l\ k = 5 -{t-ai) and (11) 
G{z) ^ {zt-H p )-\ 

H p splits the Hilbert space into a degenerate low en- 
ergy eigenspacc £_ = span{|si.Sj)|0)|V.Si, Sj G {0, 1}}, in 
which qubit k is |0), and a 5 energy eigenspacc £ + = 
span{|siSj) |1) | Vsi, Sj S {0, 1}}, in which qubit k is |1). 

First, we give the ZZXX Hamiltonian which pro- 
duces an effective a z a z interaction in the low energy 
subspace. Let Y be an arbitrary ZZXX Hamiltonian 
acting on qubits i and j and consider a perturbation 
V = Vi + V% + V3 that breaks the £_ zero eigenspacc 
degeneracy by creating an operator 0(e) close to iJ e ff in 
this space: 

Vi = [r + D(o£ + l)]®l fc -Aof<8|0}<0| fc 

V 2 = B(a x + t)®o x k (12) 

V 3 - C<j z ® |l)(l| fe . 

The term V2 above allows the mediator qubit k to un- 
dergo virtual excitations and applies an a x term to qubit 
j during transitions between the £_ and £+ subspaces. 
During excitation into £+ , the term V3 applies a a z term 
to qubit i. This perturbation is illustrated in figure 1 




FIG. 1: The ZZXX gadget used to approximate the operator 
crfa® using only a x a x and a z a z interactions. The present 
figure presents a diagrammatic representation of the Pertur- 
bation Hamiltonian V = Vi + Vi + V3 from Eq. (12) applied 
to qubits i,j and k. Not shown in the present figure is an 
overall constant energy shift of D. 

Let n± be projectors on C±] for arbitrary operator 
O we define 0± T = n±On T (0±± = n±On±) and let 
A(O) denote the lowest eigenvalue of O. One approxi- 
mates A(-fftarg) of the desired low energy effective 2-local 
Hamiltonian by a realizable 2-local physical Hamiltonian 
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H = Hp + V, where X(H) is calculated using perturba- 
tion theory. The spectrum of H is approximated by 

the projection of the self-energy operator for real- 
valued z which has the following scries expansion: 



In practice there will always be some interaction be- 
tween qubits i and j. Wc assume Y is a ZZXX Hamil- 
tonian and express the dressed Hamiltonian Y in terms 
of modified coupling coefficients. Writing the physical 
Hamiltonian: 



£__(z) = H p _ + K__+V_+G ++ (z)V + 
+ V-+G ++ (z)V + G ++ (z)V+- 



(13) 



+ o(\\v\\ 4 s 

Note that with our penalty Hamiltonian H = 0, 

and for the perturbing Hamiltonian V = V\ + V<i + V3 
only V\ is nonzero in the low energy subspace, V\ and 
V3 are nonzero in the high energy subspace, and only 
V2 induces transitions between the two subspaces. The 
non-zero projections are: 

Vi— = [Y + Aaf + D(&j + 1)]® |0> <0| fe 
V 2 _+ = B(o* + l)®|0)(l| fc 
V 2+ - = S(<7^ + l)®|l)(0| fc (14) 
V 3++ = V 3 

V+ = (Y + Ccrf+D(crJ + l))® |l)(l| fc 
The series expansion of the self-energy follows directly: 



I s * : (Y - Aa\ + D(of + 1)) d 
1 « + l) 2 ®|0)(0| fc 



yiid 



> rd 



z-S y 3 

^^(a? + lK(a| + l)®|0)(0| fe 
(o? + + 1) ® |0)(0| fe 

2 (a] + I) 3 ® 



(15) 



ADB 2 



(z-S) 

The self-energy in the low energy subspace (where 
qubit k is in state |0)) is therefore: 



£__(*) ~ Y 



A) af 



2B 2 C 
Jz~5) 2 
2B 2 4DB 2 

2B 2 C . 



(16) 



(z-sy 



(z-S) 2 1 J 

+ o (\\v\\ 4 s- 3 ) + ••• 

Y is the interaction between qubits i and j which is the 
original physical interaction dressed by the effect of vir- 
tual excitations into the high energy subspace. 



Y = Y 



B 1 



(z - S) 



Y = hiof + hj<r z + A 4 of + Aj-crJ 



JijO z o z + Kijafaj. 



(18) 



The new dressed coupling strengths are: 

2B 2 



hi ^- 

A, 



(z-S) 2 
2B 2 



(19) 



2B 2 



K, 



^t 1 ' (z-sy 

2B 2 



{z-S) 2 J ' (z-5f KlJ 
2B 2 



Kij I 1 ' (z-S) 2 J ' (z-S) 2 ' 



2B 2 



with additional couplings: 



2B 2 

:A,-1 + ^/uofof 



(z-S) 



(z-S) 



(20) 



We see that the effect of the gadget on any existing phys- 
ical interaction is to modify the coupling constants, add 
an overall shift in energy, and to add a small correction 
to the a z a z coupling which depends on the strength of 
the of term in Y . If Y is regarded as the net uncontrolled 
physical Hamiltonian coupling i and j (a source of error) 
it is only the local of field which contributes to an error 
in the a z a z coupling strength. 

We make the following choices for our gadget parame- 
ters A : B, C and D: 



A 



(21) 



2/3 



c = a -i (! 



2/3 



D 



2^1/3^2/3 



Where E is an energy scale parameter to be fixed later. 

We expand the self-energy (16) in the limit where z is 
constant (z = 0(1) < S). Writing (z - S)- 1 ~ -| + 

S ives: 



(22) 



£(0)__ = Y 

£4/3 

+ o (wvfs- 3 ) + ■ ■ ■ 

For the self-energy to become 0(e) close to Y + 
a ij a i a j ® |0)(0|fc, the error terms in (22) must be 
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bounded above by e through an appropriate choice of d. 
Define a lower bound on the spectral gap S as an inverse 
polynomial in e: 8 > Ee~ r , where E is a constant and 
integer r > 1. Now bound r by considering the (weak) 
upper bound on ||V||: 



II V1| < ||r|| + |a l3 |+4 ( 5 1 / 3 J B 2 / 3 



(23) 



E 



2/3 



E 



2/3 



Ba l a l 




The largest term in <S- 3 ||V"|| 4 is 0{E(E/S)^ 3 ), and so 
in order that <5 -3 ||V^|| 4 < e wc require r > 3. This also 
bounds the term below fourth order, i? 4 / 3 ^ 1 / 3 = O(Ee) 

and so for z <C S we obtain ||S (z) — H e s\\ = 0(e). In 

fact, S(0) = H cS + Ee(a? + 1). Now apply Theorem 

(3) from [14] and it follows that \X(H cff ) - X(H)\ = 0(e). 
It also follows from Lemma (11) of [14] that the ground 
state wavefunction of i? e ff is also close to the ground state 
of our gadget. 

The ZZXX Hamiltonian (1) allows for the direct real- 
ization of all terms in (10) except for a z a x and <j x a z 
interactions. These terms can be approximated with 
only 0(e) error using the gadget in the present section — 
thereby showing that the ZZXX Hamiltonian can effi- 
ciently approximate all terms from (10). Similarly, the 
ZX Hamiltonian allows for the direct realization of all 
terms in (10) except for a z a z and <j x <j x interactions. 
These terms will be approximated with only 0(e) er- 
ror by defining gadgets in the coming sections — showing 
that the ZX Hamiltonian can also be used to efficiently 
approximate all terms from (10). 



C. The ZZ from ZX gadget 

We approximate the operator ftijcrfa*- using the ZX 
Hamiltonian in Eq. (2) by defining a penalty Hamiltonian 
as in Eq. (11). The required perturbation is a sum of 
terms V = V\ + V 2 : 



FIG. 2: The ZZ from ZX gadget: The present figure presents 
a diagrammatic representation of the Perturbative Hamilto- 
nian V — Vi + V2 from Eq. (24) applied to qubits i,j and k. 
In addition to these terms shown in the present figure, there 
is an overall energy shift of A/2. 



The series expansion of the self-energy follows directly: 



1 s 

2 nd 

nrd 



(Y + At) 
B 2 (a z -o 



(z-S) 
B 2 



(26) 



(z-sy 



(o z -o z ]Y(a z -o z ) 



Note that in this case the desired terms appear at second 
order in the expansion, rather than at third order as was 
the case for the ZX from ZZXX gadget. The terms which 
dress the physical hamiltonian Y coupling qubits i and j 
appear at third order. The series expansion of the self- 
energy in the low energy subspace is: 



S(z)__ = (Y + At) 

2B\\-g z g z ) 
(z-S) 

+ o(\\v\\ i s- 3 ) 



(27) 



where the dressed interaction Y is defined: 
B 2 



Y = Y 



{z-S) 



: (a z -a z )Y(a z -a z 



(28) 



vi = y + A|o)(o| fc 

V 2 = B(ai-a?)(g 
The non-zero projections are: 



Vi++ 
Vi— 

V2+- 

v 2 -+ 



Y(g 
(Y- 

B(a z 



- At <E> 

<T Z 



®|1)<0|* 
® |0)(l| fc 



(24) 



(25) 



Vi does not couple the low and high energy subspaces and 
V2 couples the subspaces but is zero in each subspace. 



We assume that the physical interaction Y between i and 
j qubits is a ZX Hamiltonian and express the dressed 
Hamiltonian in terms of modified coupling constants. 
Writing the physical Hamiltonian: 



Y 



hiU z + hj<r z 

+ J l3 <T Z <T X 



K v] a x a z 



We obtain modified coupling strengths: 



hi 
h ; 



hi 
hi 



2B 2 (h l - h 3 
(z - S) 2 

2B 2 (h ] - hi 
(z - 6) 2 



(29) 



(30) 
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In this case only the local Z field strengths are modified. 
We choose values for the perturbation interaction 



and A 



fin and ex- 



strengths as follows: B — y — ^ — HtJ 
pand the self-energy in the limit where z is constant 
(* = 0(1) « 



S(0)__=F + Ai^l 
+ 0(||V r || 4 (5- 3 ). 



(31) 



We again choose 8 to be an inverse power in a small 
parameter e so that <5 > Ee~ s , and again use the (weak) 
upper bound on \ \V\\: 



< 



\\Y\\+Pi: 



(32) 



fj8 1 , and so in order 



The largest term in ||V|| 4 (5- 3 is 0^5 
that ||V|| 4 <5~ 3 < e we require r > 1. 

Using the gadget defined in the present section, the ZX 
Hamiltonian can now be used to efficiently approximate 
all terms in (10) except for a x a x interactions. These 
interactions can also be approximated with only 0(e) 
error by defining an additional gadget in the next section. 



D. The XX from ZX gadget 

An a x a x coupling may be produced from the a z a x 
coupling as follows. We define a penalty Hamiltonian 
and corresponding Green's function: 



H p 

G ++ 



z-S' 



s\-)(- 



-)(- 



(33) 



This penalty Hamiltonian splits the Hilbert space into a 
low energy subspace in which the ancilla qubit k is in 
state |+) = (|0) + |l))/\/2 and a high energy subspace in 



|1»/V2- 

Vl + V 2 , 



(34) 



(35) 



which the ancilla qubit k is in state |— ) = (|0) 
The perturbation is a sum of two terms V - 
where Vi and V 2 are given by: 

Vx = Y®l k + A\+)(+\ k 
V 2 = B(af - a x )al 

The non-zero projections are: 

Vi ++ = l-X-N-X-l* 

= Y®\-){-\ k 

Vi— = Y®\+){+\ k +A\+)(+\ k 

V 2+ - = \-)(-\V 2 \+)(+\k 

= B(a x -a x )\-)(+\k 

V2-+ = \+)(+\V2\-)(-\k 

= B (af-a?)\+)(-\ k . 

Once more we see that the perturbation V\ does not cou- 
ple the subspaces, whereas V 2 couples the subspaces but 




FIG. 3: The XX from ZX gadget: The present figure 
presents a diagrammatic representation of the Perturbative 
Hamiltonian V = V\ + Vi from Eq. (34) applied to qubits 
i,j and k. In addition to the terms shown in the present fig- 
ure, there is an overall energy shift of A/2. The penalty term 
applied to qubit k is the a x basis. 



is zero in each subspace. This perturbation is illustrated 
in Figure 3. 

The series expansion of the self-energy follows: 



I s * : (Y + Al)®\+)(+\ k 



(36) 



3 rd 



B 2 {o x -a x f 

{ 1 ]> ol+X+U 



(z-S) 
B 2 



(z-S) 



(a x -a x )Y(a x -a x ) 



2 V J 



Again we see that the desired term appears at second 
order, while the third order term is due to the dressing 
of the physical interaction Y between qubits i and j. In 
the low energy subspace the series expansion of the self- 
energy to third order is: 



S(z)__ = (Y + At) 

2B 2 (l-a x a x ) 
+ {z~S) 

+ o(\\v\\ 4 s- 3 ) 

where the dressed interaction Y is defined: 
B 2 



(37) 



Y = Y 



{z-S)- 



(of - a x )Y(af - of). (38) 



Once more we assume the physical Hamiltonian Y is a ZX 
Hamiltonian 29 and we describe the effects of dressing to 
low order in terms of the new dressed coupling strengths: 



2B 2 {A l -A J ) 

(z - 5) 2 
2B 2 (A J -A t ) 

(z - 5) 2 ' 



(39) 



and in this case only the local X field strengths are mod- 
ified. 

Choosing values for our gadget parameters A — 7^ 
and B = J 2^ and expanding the self-energy in the 
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limit where z is constant (z — 0(1) <C S) gives: 

S(0)__ = Y®\+)(+\ k (40) 
+ JijafaJ ®|+)(+U 
+ 0(||F|| 4 <T 3 ) 

As before, this self-energy may be made 0(e) close to the 
target Hamiltonian by a bound 5 > Ee^ 1 . 

b. Summary The proof of Theorem (1) follows from 
the simultaneous application of the ZZXX gadget illus- 
trated in Fig. 1 to realize all a z a x terms in the target 
Hamiltonian using a ZZXX Hamiltonian. Similarly, ap- 
plication of the two gadgets illustrated in Fig. 2 and 3 
to realize a x a x and <r z a z terms in the target Hamilto- 
nian proves the first part of Theorem (2). Our result is 
based on Theorem (3) from [14] which allowed us to ap- 
proximate (with 0(e) error) all the Hamiltonian terms 
from Eq. (3) using either the ZZXX or ZX Hamiltonians. 
It also follows from Lemma (11) of [14] that the ground 
state wavefunction of i? e ff is also close to the ground state 
of our gadget. So to complete our proof, it is enough to 
show that each gadget satisfies the criteria given in The- 



orem (3) from [14]. 



II. CONCLUSION 

The objective of this work was to provide simple model 
Hamiltonians that are of practical interest for experi- 
mentalists working towards the realization of a univer- 
sal adiabatic quantum computer. Accomplishing such as 
task also enabled us to find the simplest known QMA- 
complete 2-local Hamiltonians. The a x a x coupler is re- 
alizable using systems including capacitive coupling of 
flux qubits [22] and spin models implemented with po- 
lar molecules [23]. In addition, a a z a x coupler for flux 
qubits is given in [24] . The ZX and ZZXX Hamiltonians 
enable gate model [25], autonomous [26], measurement- 
based [27] and universal adiabatic quantum computa- 
tion [10, 14, 15], and may also be useful for quan- 
tum annealing [28]. For these reasons, the reported 
Hamiltonians are of interest to those concerned with the 
practical construction of a universal adiabatic quantum 
computer [30]. 
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